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Nonlinear Systems 


x=Ax+Bu —> x= f(x,u) (1) 
Analysis: We use the shorthand notation x = f(x) 
for 4 x(t) = f(x(t)). 
x=f(x) f:R’>R" time-invariant (autonomous) 
x=f(t,x) f:IRx R"— R" time-varying (non-autonomous) 


Design: 


x = f(x,u) u to be designed as a function of x. 
Equilibria 
x = x* is an equilibrium for x = f(x) if f(x*) =0. 


Example: Linear system x = Ax. 
If A is nonsingular, x* = 0 is the unique equilibrium. 


If A is singular, the nullspace defines a continuum of equilibria. 


Example: Logistic growth model in population dynamics 


z= f(x)=r(1-2)x, r>0, K>0 (2) 
= 


growth rate 


x > 0 denotes the population and K is called the carrying capacity. 


For systems with a scalar state variable x € R, stability can be deter- 
mined from the sign of f(x) around the equilibrium. In this example 
f(x) > 0 for x € (0,K), and f(x) < 0 for x > K; therefore 


x=0 unstable equilibrium 


x=K_ asymptotically stable. 
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Linearization 


Local stability properties of x* can be determined by linearizing the 
vector field f(x) at x*: 


a 
FX 2 SF le of ¥ + higher order terms (3) 
—— ox x=x* 
=A 


Thus, the linearized model is: 


% = AX. (4) 


If RA;(A) < 0 for each eigenvalue A; of A, then x* is asymp. stable. 


If RA;(A) > 0 for some eigenvalue A; of A, then x* is unstable. 


Example: Logistic growth model above: 


f(x) 


\ 
f'(0) >0  f'(K) <0 
unstable stable 


Caveats: 


1. Only local properties can be determined from the linearization. 
Example: The logistic growth model linearized at x = 0 (x = rx) 
would incorrectly predict unbounded growth of x(t). In reality, 
x(t) — K. 

2. If RA;\(A) < 0 with equality for some i, then linearization is 
inconclusive as a stability test. Higher order terms determine 
stability. 

Example: fo=s vs. f(x) =—-x3 


f'(0) = 0 in each case, but one is stable and the other is unstable. 


2 
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Second order example: Pendulum 


(m6 = —kl6 — mg sin@ (5) 
Define x = ; } State space: S! x R. 
Xy = X2 
k os (6) 
x2 = ——X2—- 7 SIN X, 
m £ 
Equilibria: (0,0) and (71,0) 
0 
g ‘ (stable) at x; = 0 
of a 0 1 —F 79 
ox —Scosx; —— | 0 1 
E : (unstable) at x; = 71 
£  ¢ 


5 SSS = SS iam = Z ome Figure 1: Phase portrait of the pendu- 
. lum for the undamped case k = 0. 


Essentially Nonlinear Phenomena 


1. Finite Escape Time 


Example: * = x? 
a 
ee R= 1 
1 1. ‘ 
x(t)  x(0) (7) 
+ 4oo—= 
x0) 


For linear systems, x(f) —> co cannot happen in finite time. 


jn} -------- 55S TSS. 


z 
= 
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2. Multiple Isolated Equilibria 
Linear systems: either unique equilibrium or a continuum 
Pendulum: two isolated equilibria (one stable, one unstable) 
“Multi-stable” systems: two or more stable equilibria 
Example: bistable switch 


xy = —ax, + x2 x1: concentration of protein 


xt (8) 


Xo = — bx. Xx: concentration of mRNA 


a >0,b > Oare constants. State space: R>o x Rso. 


This model describes a positive feedback where the protein en- 
x 


z+ 
1+x7 


Single equilibrium at the origin when ab > 0.5. If ab < 0.5, the line 


coded by a gene stimulates more transcription via the term 


where x; = 0 intersects the sigmoidal curve where x2 = 0 at two 
other points, giving rise to a total of three equilibria: 


x2 A 


> XxX 


stable unstable stable 
(gene off) (saddle point) (gene on) 


